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Abstract We study the mechanism of formation of Bose glass (BG) phases in the spin-1 Bose Hubbard 
model when diagonal disorder is introduced. To this aim, we analyze first the phase diagram in the zero- 
hopping limit, there disorder induces superposition between Mott insulator (MI) phases with different filling 
numbers. Then BG appears as a compressible but still insulating phase. The phase diagram for finite hopping 
is also calculated with the Gutzwiller approximation. The bosons' spin degree of freedom introduces another 
scattering channel in the two-body interaction modifying the stability of MI regions with respect to the action 
of disorder. This leads to some peculiar phenomena such as the creation of BG of singlets, for very strong 
spin correlation, or the disappearance of BG phase in some particular cases where fluctuations are not able to 
mix different MI regions. 

PACS numbers: 03.75.Mn,64.60.Cn,67.85.-d 



1 Introduction 

Ultracold atomic gases can be described by interacting atoms with an internal spin, corresponding to a low- 
energy hyperfine level F. If the spin orientation is fixed by an external magnetic field, as it happens when the 
gas is confined in magnetic trap, a scalar model is sufficient to describe the system. Conversely, if the spin 
orientation is not externally constrained, as in the case of optical trapping, the spinor character of the gas has 
to be taken into account. For bosons trapped in a deep optical lattice potential the system is well described 
by the spinor Bose-Hubbard (BH) model 1 1]. Bosonic interactions, treated as two-body contact collisions, are 
sensitive to the spin degree of freedom and contribute to the orderings at zero temperature. 

As in the scalar case [2], the competition between hopping and interactions leads to a quantum phase 
transition between spinor superfluid (SF) condensate and a Mott insulator (MI) state IHIllH]- Spin correlation 
introduces magnetic ordering which contribute to the stability of one phase with respect to the other. Moreover, 
the presence of spin scattering channels influences also the stability of the MI phases in the presence of 
different types of disorder. 

Disorder plays an essential role in condensed matter physics and it has been shown to be an essential 
ingredient for studies of conductivity, transport, high-Tc superconductivity, neural networks or quantum chaos 
to mention few examples (see the reviewlSl] and references therein). Disorder can be produced in ultracold 
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atoms in a controlled and reproducible way. Standard methods to achieve such a controlled disorder are the 
use of speckle pattems SIH] which can be added to the confining potential, or optical superlattices created by 
the simultaneous presence of optical lattices of incommensurate frequencies lisilol flOll . Other methods include 
using an admixture of different atomic species randomly trapped in sites distributed across the sample and 
acting as impurities 1 1 1, 12], or the use of inhomogeneous magnetic fields which modify randomly, close to a 
Feshbach resonance, the scattering length of atoms in the sample depending on their spatial position lfl3L[l4ll . 

Recently, the phase diagram of the spin-1 BH model in two dimensions (2D) in the presence of disorder has 
been studied with a Gutzwiller mean field approximation! 15]. As in the scalar case [2], a gapless Bose-glass 
(EG) insulator phase appears, characterised by finite compressibility and exponentially decaying superfluid 
correlations in space, but, because of the spin interaction, the phase diagram changes considerably. 

In this paper we focus on the spin-1 BH model in the presence of diagonal disorder and analyse the role 
of spin correlations in the formation of the EG phase. We study in detail the zero-hopping limit (atomic case). 
We provide, in this simple case, the phase diagram corresponding to different types of disorder and show 
how spin correlations can prevent the formation of the BG between some MI regions. This fully analytical 
approach allows to easily visualise the mechanism of formation of the BG in terms of superpositions between 
MI phases with different filling factors. We also provide, in some cases, the complete phase diagram for 
finite tunnelling using a numerical mean field Gutzwiller approximation. In particular, we show the case 
of large spin interactions where a BG of singlets emerges. We analyse the case in which the disorder is 
directly introduced in the two-body scattering lengths, associated with total spin of scattering particles s = 
and = 2, instead of in the Hamiltonian parameters directly. This represents a more realistic scenario since 
small fluctuations in the scattering lengths can be introduced by optical Feshbach resonances. We observe the 
absence of BG phase for ferromagnetic spin interactions and disorder in the s = 2 scattering length, confirming 
that, in this limit, the spinor model is equivalent to the scalar case. 

The paper is organised as follows: in Sec.|2]we introduce the BH model. Then we proceed to compare 
the effect of diagonal disorder in the local potential or in the interaction. In particular, we compare the cases 
in which disorder is in the local parameters of the Hamiltonian ( Sec. O with the case in which disorder is 
introduced directly in the two-body scattering lengths corresponding to the collision channel with total spin 
s — and s — 2 (Sec.|4l). Finally, in Sec.|5]we present our conclusions. 



2 Model 

Low energy spin-1 bosons loaded in optical lattices, sufficiently deep so that only the lowest energy band is 
relevant, can be described by the spinor BH model. The corresponding Hamiltonian is UJ: 

where (ij) indicates that the sum is restricted to nearest neighbours in the lattice and d]^ (a,(j) denotes the 
creation (annihilation) operator of a boson in the lowest Bloch band localised on site / with spin component 
a = 0,±l. 

The first term in ([T} represents the kinetic energy and describes spin independent hopping between nearest- 
neighbour sites with tunnelling amplitude f . The second and third term account for spin independent and spin 
dependent on site interactions, respectively. These energies at site / are defined as Uq,2 = cq ^ / drw'^{r — r,) 
with Co = 47r7i^(ao +2fl2)/ (3m) and C2 = Anfi^{a2 — an) I (3m), where as with 5 = 0, 2 is the s-wave scattering 
length corresponding to the channel with total spin S Ilia. [ITIi and w{r — r,) is the Wannier function of the 
lowest band at site While the second term of ^ is spin independent and equivalent to the interaction energy 
for scalar bosons, the third term represents the energy associated with spin configurations within lattice sites 
with 

S/= L ^Gi^co'ac'i, (2) 

(J(J'=0,±1 

being the spin operator at site / and F the traceless spin-1 matrices. S's components obey standard angular 
momentum commutation relations = iSijkSk- The spin-interaction term favours a configuration with 

total magnetisation zero for U2 > 0, denoted as polar and sometimes antiferromagnetic. The ferromagnetic 
configuration where spins add to a maximal possible value corresponds to U2 < llT6l[l7l[l8ll. In the grand 



^hi{hi -~l) + ^{Sf - 2hi) - nhi 



(1) 
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canonical approach the total number of particles is controlled by the last term of (HJ where /i is the chemical 
potential and 

CJ=0,±1 

is the total number of bosons on site i. Hamiltonian ([ij can be straightforwardly derived from the micro- 
scopical description of bosonic atoms, with a hyperfine spin F = I, loaded in a deep optical lattice and 
considering the two-body short range (s-wave) collisions. More details about the derivation can be found in 

lUIillnl El Si- 
Notice also that, since the orbital part of the wave function in one lattice site is the product of Wannier 
functions for all the atoms, it is symmetric under permutation of any two atoms. Therefore, the spin part of 
the wavefunction should also be symmetric due to Bose statistics. This imposes .s, +«, to be even 1I21I1 . being 
Si and «, the quantum numbers labelling the eigenvalues of S, and «,, respectively. In alkaline atoms usually, 
the scattering lengths are similar, — 02, and the symmetry of the Hamiltonian becomes SU(3) instead of 
SU(2). That implies |C/o| > |C/2l- 

As in the scalar case, the spinor BH system exhibits a quantum phase transition between superfluid and 
insulating states id Si- In the insulating states, fluctuations in the atom number per site are suppressed and 
virtual tunnelling gives rise to effective spin exchange interactions that determine a rich phase diagram in 
which different insulating phases differ by their spin correlations. 



2.1 Atomic limit (f = 0) 

Hamiltonian ([T]) becomes diagonal in the limit of vanishing hopping. In this limit it reduces to the sum of 
local terms H — > with 

^0 = Yn(n-l) + ^(s2-2n)-M«, (4) 

(site indexes are omitted if the system is homogeneous). Eigenstates of Hq are denoted with |n,5,m), where 
the quantum numbers refer to the three commuting local observables 

n|n,i,m) = n|n,i,m) , 

|n,i,m) = s{s+ 1) \n,s,m) , 
S, |n,5,m) = m\n,s,m) , 

with energies //q jn, 5, m) = Eo{n,s)\n,s,in) given by 

Eoin,s) = -iin + - 1) + y {s{s+\) - 2n) . (5) 

From Eq.(|5]l, one can easily deduce the structure of the ground state of the insulator phases in the limit 
t — the atomic limit without disorder, only MI phases exist. SF corresponds to the points separating MI 
intervals with different fillings (as can be inferred from f — > limit). This means that, at fixed /i, the ground 
state has an integer filling n. The boundaries corresponding to the degeneracy points between two fillings n\ 
and n2 satisfy the condition i?o(«ij'Si) = £^o(«2)'*2)- Any MI region with filling n and fixed U2 is defined for 
li-{n) < /i < where /i±(n) are the boundaries of the region. 

For antiferromagnetic interactions, U2 > 0, the minimum energy Eq'" is attained with minimum s, its 
specific value depending of the number of atoms per site. Thus, for the even filling factor, the minimum spi n 
is zero and the state is described as [0,0, «) with n even. This state is known as spin singlet insulator ll22ll . 
If the atom number per site is odd, then the minimum spin per site is one and the state reads \l,m,n) in the 
absence of disorder. 

Let us calculate explicitly jj.±. For C/2 > (antiferromagnetic case) we must distinguish between odd and 
even occupation lobes: 

(1) n odd: the boundaries /i±(«) are obtained by imposing EQ{n, 1) = Eo{n± 1,0) 

11^ = («- l)[/o, 

^+ = nUo-2U2. (6) 



Odd lobes exist only for jU_ < ^+ that is, for U2/U0 < 0.5. 
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Fig. 1 Phase diagram of the spinor F = 1 BH model in the limit t = 0. Each region corresponds to a MI phase with a dif- 
ferent occupation number. For U2/UQ > the system has zero magnetization (antiferro) and for U2/U0 < the system is in a 
ferromagnetic phase 



For U2 <0 (ferromagnetic case) the maximum value of the spin is given hy s = n and the condition for 
the boundaries £0(11, n) = £0 (« T 1 j «) reads 



The atomic-limit phase diagram is depicted in Fig. [T] the MI intervals obtained fixing U2 correspond to 
the basis of the MI lobes in the t/Uo - Ix/Uq plane lIlSll . In the antiferromagnetic region, one can see that the 
formation of singlet stabilises the even MI lobes while the odd lobes shrink. For U2 < 0, U2, eq. ([T| reduces 
to the scalar Hamiltonian with Uq + U2 put in place of Uq. Notice that for U2 < ~Uo the spectrum (|5]l is not 
bounded from below and the model becomes instable. 



3 Diagonal disorder in the Hamiltonian's parameters 

Starting from this scenario, we proceed to analyse the stability of the MI phase in the presence of disorder. 
We start by analysing local fluctuations e, added either to the homogeneous chemical potential ji, or on 
the interaction potential Uq or 1/2- We consider e, to be a random variable defined for every site i with a 



(2) n even: for U2/U0 < 0.5, H±{n) are obtained by imposing EQ{n,0) = Eo{n±l,l) 

II- = {n-l)Uo-2U2, 
jx+ = nUo, 

while, for U2/Uo> 0.5 the condition £0(0, n) = Eq{0, n =f 2) leads to 



(7) 




(8) 



= (n-l)(f/o + C/2), 
^i+=niUo + U2). 



(9) 
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given probability distribution p{£). Here we consider a bounded probability distribution that is —A < £, < 
A. Disorder mixes together different MI regions so that degeneracy between different fillings appears. So, 
between MI intervals, regions can appear where filling is not defined and correspondingly a BG phase appears. 
Being the probability distribution bounded, the BG phase appears only around the original boundaries and its 
extension depends on A . On the other hand, deep inside the MI, small fluctuations of the disordered parameter 
are not able to mix different fillings and the state remains stable in the corresponding MI phase. 

The first case we study is the disorder in the chemical potential. This type of disorder can be produced 
introducing some random inhomogeneities in the local potential fij = {1 + Ej. Fixing U2, disorder in /i^ cor- 
responds to horizontal fluctuations of maximum amplitude A in the atomic-limit phase diagram. The new 
boundaries for the MI phases are given by replacing /i± — > /i± =F Zi as can be observed in the top panel of Fig. 
[21 We notice that BG regions always appear between MI lobes with width (in /i) of the order of 2 A . Since odd 
lobes shrunk by U2, they are more unstable and disorder can make them disappear. 

When U2/U0 > 0.5 disorder mixes together only even occupations so BG is formed by singlets llSll . The 
finite-hopping phase diagram for this case is displayed in the bottom panel of Fig. |2] where we have used 
the condensate fraction as an order parameter to separate SF from BG phase. Our calculations are done 
within the mean field Gutzwiller approach, which is able to take into account inhomogeneties caused by 
disorder!. 15.1. The MI phase is characterised by vanishing fluctuations in the density (or zero compressibility) 
while the BG corresponds to finite density fluctuations and zero condensate fraction. In this plot one can see 
the disappearance of the odd occupation MI lobes. Numerical calculation shows that (S^) = meaning that 
also the BG phase is formed by singlets 1 15]. Note that, even if we found it adding disorder in the chemical 
potential, singlet BG appears also for disorder in Uq and U2, the only necessary condition being U2/U0 > 0. 

Disorder in U2 corresponds to fluctuations along the vertical direction of the phase diagram in Fig. [T] 
The resulting phase diagram in the atomic limit is shown in the top panel of Fig. [3] The new boundaries are 
obtained changing jj.± in Eqs. replacing C/2 — > C^2 for U2/U0 < 0, and U2 ^ U2±A for U2/U0 > . 
One consequence is that, for U2/UQ < 0.5 no BG appears between even and odd lobes with higher occupation 
(that is between lobes with occupation « = 2m and n = 2m + 1) where the boundaries are vertical. This is 
true for \U2\ > A otherwise disorder can mix together ferromagnetic and anti ferromagnetic regions, and BG 
appears between all the lobes. Note that, in the antiferromagnetic limit, the shrinking of the MI lobes is almost 
independent on the occupation and depends only on U2- In th ferromagnetic case, the MI lobes become more 
unstable increasing n disappearing for for n> {U0 + U2 + A)/ {2A ) . 

So far we have analysed the effects of fluctuations along vertical and horizontal directions in the phase 
plane. More generally, a diagonal disorder can be interpreted in terms of fluctuations of a certain amplitude A 
along some direction, possibly depending on the position (/i^jt/j) in atomic phase diagram, in which is 
centred. Let us denote the direction as a straight line 

^=a-^+b, (10) 

In general, if a MI boundary lies along the direction of fluctuations defined by a specific type of disorder, 
no BG would appear for sufficiently small disorder strength. 

In the case of disorder in Uq, the fluctuations occur along the direction defined by the parameters in (llOb 



a = U^/jX^ and = with amplitude A=2A^ {ji^)^ + {U^Y/{Uq - A^). The new boundaries are obtained 
modifying pi± by the replacements Uq ^ Uq^ A. This causes a progressive shrinking of all the lobes which 
disappear for high densities. It is worth noting that odd lobes disappear independently on U2 for n > ^^^^ • 
The resulting atomic phase diagram is shown in the bottom panel of Fig. [3] 



4 Disorder in the scattering lengths 

As pointed out in Sec.|2] Uq and U2 depend on the scattering lengths aq and 02 via the relations 

Uq = ao + 2a2, 

[/2 = ao-a2, (11) 

having introduced the renormalized parameters Us = a^^K J drw'^{r — r,) /(3m). Both Oq and a2 can fluctuate 
locally in the presence of optical Feshbach resonances, so disorder can be introduced in these variables and 
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Fig. 2 (Color online) Top panel: phase diagram in the limit t = and disorder in /x with A /Uo = 0.4. Bottom panel: condensate 
fraction with disorder in fi for Ui/Uq = 0.7 and A /Uo = 0.4. The solid lines correspond to the boundaries of the MI lobes. The 
region outside the MI phase and with vanishing condensate fraction corresponds to the BG. In this case BG is formed by singlets. 



the new scattering lengths become site-dependent a/ = + £/. This represents a more realistic scenario than 
considering disorder in Uq or U2 alone. 

To find the phase diagram in the atomic limit, we have to modify again the boundaries in (|6]|9]l as we 
did in the previous section in the cases of disorder in Uq and U2, remembering that, this time, both of them 
fluctuate at the same time. So, in the case of disorder in a2, one has to do the replacement Uq ^ UqT 24 
and U2 ~^ U2T ^ for U2/UQ > and U2 — > i72 ± ^ for U2/U0 < 0. The fluctuations occur, with amplitude 

A=2A ^J{2|lY + {Uo - 2U^f/{UQ - AA^), along the direction ^ with a = {U^ + Uq) /jM) and b = 0.5 
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Fig. 3 Top panel: phase diagram in tiie limit t = and disorder in U2 with A /Uq = 0.1. Bottom panel: phase diagram in the limit 
t = and disorder in Uo with A /Uo = 0. 12. 



which is a family of straight lines passing through {fi/Uo = O^Uj/Uq = 0.5). This means that the boundary 
between first and second lobe belongs to this family and no BG is expected between these two MI regions, as 
can be observed in both the atomic-limit (top panel) and complete phase (bottom panel) diagram in Fig.|4] 

If the disorder is set in Oq, substitutions in Eqs. (|6]|9]| are Uq Uq^ A and Uj^ U2±A for U2/UQ > 
and Us UsT^ (■* = 0, 1) for Uj/Uq < 0. The direction of the fluctuations is given by a = (C/j + ^0) /Mo 

and = -1 and amplitude A=2A ^ {pL^f + {UQ + U^f /(Uo - A'^). In this case, all the boundaries in the 

ferromagnetic limits lie on fluctuation directions so no BG appears in the ferromagnetic regime as is shown in 
the top panel of Fig. |5j This result marks a distinction between the scalar case and the spinor one with ferro- 
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Fig. 4 (Color online) Disorder in a2 with A /Ug = 0.04. Top panel: atomic limit. Bottom panel: finite hopping phase diagram 
with (72 /{/o =0.1. 

magnetic spin correlations where disorder in the s = scattering channel only is not enough to produce BG. 
The finite-hopping phase diagram for the antiferromagnetic regime where MI lobes are always surrounded by 
BG is shown in the lower panel of Fig|5] 



5 Conclusions 



We have analysed the spin-1 Bose Hubbard model with different types of diagonal disorder, focusing on the 
atomic limit to illustrate how Bose Glass phase emerges in between Mott lobes. In this limit, disorder mixes 
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Fig. 5 (Color online) Disorder in Oq with A /Ug = 0.04. Top panel: atomic limit. Bottom panel: finite hopping phase diagram 
with (72 /[/o =0.1. 



MI phases with different occupation numbers, producing regions of BG between MI intervals. The study of 
the atomic limit gives useful information also about the finite hopping case, providing the structure of the 
phase diagram near the basis of the MI lobes. To illustrate the power of this approach, we have also shown in 
some cases the complete phase diagram, calculated by Gutzwiller approximation. We first analysed disorder, 
in either the chemical potential or the local interaction Uq and U2, explaining how to construct the phase 
diagram in the atomic limit. Then, we used these results to study the more realistic case of disorder in one 
of the two scattering lengths aq or fl2 corresponding to different scattering channels. While in the scalar case 
BG always appears between MI regions, as soon as a small disorder is introduced, the spinor character can 
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stabilise the MI phase for some densities or spin interactions inhibiting the BG creation near its boundary. 
That happens in the case of disorder in U2 between lobes with n = 2m and n = 2m + 1, for disorder in ao in 
the ferromagnetic limit corresponding to f/i/f^o < as well as for disorder in 02 between the first and second 
lobe. The creation of singlets in the even occupation MI enhances the stability of this phase, reducing the 
odd occupation MI lobes. As a consequence, the odd lobes are also less stable under the effect of disorder 
In the extreme case in which odd lobes disappear also without disorder, the BG assumes a spin structure 
of singlet. As a future perspective, we would like to extend the approach we have illustrated in this paper 
to better enter into the spin properties of BG phases as well as the MI in the regimes where ferromagnetic 
and antiferromagnetic orders are mixed by disorder without destroying the MI phase. In this limit some new 
interesting phases could appear where the glassy character is not embedded in the density but in the spin 
degrees of freedom. 
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